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Elastic backscattering of electrons moving along the helical edge is prohibited by time-reversal 
symmetry (TRS). We demonstrate, however, that an ensemble of magnetic impurities may cause 
TRS-preserving quasi-elastic backscattering, resulting in interference effects in the conductance. The 
characteristic energy transferred in a backscattering event is suppressed due to the RKKY interaction 
of localized spins (the suppression is exponential in the total number of magnetic impurities). We 
predict the statistics of conductance fluctuations to differ from those in the conventional case of a 
one-dimensional system with quenched disorder. 



The constitutive porperty of a two-dimensional topo- 
logical insulator is the presence of helical edge states at 
its bounaries. The electron states propagating in oppo- 
site directions along an edge form Kramers doublets and 
are protected against elastic backsacattering by the time- 
reversal symmetry (TRS). The fundamental consequence 
is the universality of the zero-temperature conductance 
of a topological insulator. Experimental demonstration 
of such universality is viewed as the confirmation of the 
existence of topological insulators. The existing experi- 
ments clearly distinguish the topological insulators from 
highly-resistive "conventional" onea^. However, the mea- 
sured conductance approaches the universal value only in 
very short (less than l/j,m long) samples. Conductance of 
longer samples typically is lower, indicating the presence 
of electron backscattering. 

Mechanisms of the backsacttering are a matter of ongo- 
ing debate. Current proposals utilize the Coulomb inter- 
action between the electrons of helical edg^^^ or scatter- 
ing of an electron off a localized magnetic impuritji^ls! 
In all theoretical models considered so far the electron 
backscattering is either deeply inelastic or, as in the case 
of a magnetic impurity at temperatures exceeding the 
Kondo scale, quasi-elastic and incoherent due to the flips 
of the impurity spin'^. In either case the bacskcattering, 
while suppressing the value of conductance G at finite 
temperatures does not lead to mesoscopic fluctuations in 
G, in contrast with the low-temperature electron trans- 
port in conventional low-dimensional conductors. 

In this work we show that interference effects in the 
conductance of a helical edge contaminated by an en- 
semble of magnetic impurities are actually possible. Un- 
like a single impurity, spins in an ensemble are coupled 
by a long-range RKKY interaction. The latter prevents 
individual spins from flipping in the course of electron 
backscattering. That allows for a coherent superposition 
of the electron de Broglie waves reflected by different im- 
purities. Remarkably, this picture does not contradict 
TRS: the electron reflection remains inelastic. However, 
it is associated with a collective flip of a block of spins. 



Therefore its amplitude is sensitive to the spatial struc- 
ture of the electron wave. Large number of spins in the 
block also ensures a parametrically small (exponential 
in the number of spins) energy transfer in the scatter- 
ing event. This energy sets a mild lower limit for the 
temperature at which the considered quantum-coherent 
phenomenon can be observed. 

Interference in electron reflection manifests itself in 
mesoscopic fluctuations of conductance G, seen, for ex- 
ample when the chemical potential of electrons is contin- 
uously changed by a gate voltage Vg. In conventional one- 
dimensional wires such fluctuations are due to the vary- 
ing interference conditions for the potential scattering of 
electrons. At sufficiently low temperatures such conduc- 
tance fluctuations are pronouncedly non-Gaussian. In 
the regime of weak backscattering they obey the Rayleigh 
statistics. In contrast, we find that the variations of G 
with Vg caused by an ensemble of magnetic impurities 
tend to be close to gaussian except for temperature close 
to the spin-glass crossover. 

We consider a topological insulator with a simple he- 
lical edg^^l such that the spin s of an electron occu- 
pying an edge state has a conserved component Sz in 
some fixed direction z. A magnetic impurity in a vicin- 
ity of the edge will experience two important interac- 
tions: a local single-ion anisotropy induced by the bulk 
spin-orbit coupling and a local exchange coupling to the 
electrons of the edge. (The direct exchange between the 
impurity spins is negligible at low impurity density.) We 
assume the anisotropy to be of easy-axis type with some 
anisotropy constant K > Q. The exchange coupling of the 
edge-state electrons to the impurity spin will generally be 
anisotropic and depend on the position of the impurity. 
The effective low-energy Hamiltonian describing the he- 
lical edge with TV magnetic impurities is 
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Here the two-component spinor fields ip{x) represent 
the smooth (on the scale provided by the Fermi wave 
length 2n/kp) envelope of the electron operators, t{x) = 
{x)e~'^^'''^'^Te^^''^'^'il}{x) is the electron spin density 
operator where r is the spin vector composed of the three 
Pauli matrices; v is the electron velocity, S' is the I'th 
impurity spin. We will see that the interference effects in 
G appear if the spin anisotropy axis n is different from 
z and the impurities have spin S > 1. We assume the 
magnetic impurities to be distributed randomly along the 
sample length L and at random distances from the edge, 
resulting in random positions Xi and coupling constants 
k\. In general, the exchange tensors k\^^ and the ten- 
sor of single-ion anisotropy should be considered as run- 
ning coupling constants, depending on the choice of the 
bandwidth cutoff. The renormalization of the anisotropy 
is not infrared-divergent. Therefore, assuming that the 
bare constants are small, one adds an ultraviolet cor- 
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rection to Eq. (|T|) of the form 
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A, where A is the bandwidth cutoff 



with \5K^\ 

scale set by the insulator band gap. Renormalization of 
the tensors is infrared-divergent, but remains small 
at energies exceeding the Kondo scale Tr-. The latter is 
exponentially small in 1/|k*| or even a higher power of 
that parameter, due to the effect of single-ion anisotropy. 
We note that the itinerant electrons facilitate the RKKY 
interaction between the impurities, which is of the order 
of for two impurities at distance x from each 

other. We may ignore the renormalization of k as long 
as the RKKY exchange is larger than T^- 

We start the analysis of the model with considering two 
impurities at distance x from each other. Treating the 
exchange tensors as perturbation theory parameters 
we find the leading-order RKKY interaction 
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Here Ldi,c{x) is the orthogonal matrix of counterclock- 
wise rotation through angle 2kpx about the z axis and 
Pcd = ^cd ~ ^cz^dz is the matrix of orthogonal projec- 
tion onto the xy plane. Assuming that Hrkky, Eq. Q, 
and T-La, Eq. ([s]), are small as compared to the easy- 
axis anisotropy K, we may apply secular perturbation 
theory to determine the low-energy spectrum of the two- 
spin system. In the zeroth-order perturbation theory the 
ground state of the two-spin system is four-fold degener- 
ate with the corresponding eigenspace spanned by four 
vectors, | ± 5)1 (g) | ± 8)2^ where \s)i denotes an eigenstate 
of S* -n with an eigenvalue s G {—5, — S*-!-!, . . . , S*— 1, S} 
and the subscript i = 1,2 labels the Hilbert space at- 
tached to the ith spin. The secular matrix of perturba- 
tion Eq. Q is conveniently written in terms of operators 
Si acting on the i-th spin only: 



(5) 



In this notation, the secular matrix takes form of Ising 
Hamiltonian, 



%! = —6E cos{2kpx)siS2 ; 
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here 5E = {hv / ATr\x\)naK^i,PbcK'^d''^d,- Depending on the 
sign of cos{2kFx) the ground state of the Hamiltonian 
(p^ is one of the two doublets, \S)i(S)\S)2, \- S)i(S)\- 8)2 
ot\S)i(^\-S)2, \-S)i^\S)2- The ground-state doublet is 
separated by energy ~ \SE\ from the excited level, which 
is also a doublet. Note that for all S" > 1 the contribution 
of the perturbation Eq. (|3| to the secular matrix is purely 
diagonal and has no effect on the splitting of the ground 
level. 

In higher-order perturbation theory further splitting 
of the two doubly degenerate energy levels occurs. The 
dominant effect here is due to the perturbation ([s]). In- 
deed, it has non-vanishing matrix elements for transi- 
tions in which the projection of one impurity spin is in- 
creased (or decreased) by 1 or 2, for example {{S — 2|(8) 
{S\)nA{\S)(»\S)) ^ 0. For 5 > 1 any such transition 
takes a vacuum state to a virtual state having the energy 
of the order K. At least [S— 1/2] such consecutive transi- 
tions (here the symbol [. . . ] stands for the integer part of 
S) are needed in order to flip one impurity spin from —S 
to S, that is to bring the system from a vacuum state to 
a state in the excited doublet with energy \SE\. Taking 
such processes into account amounts to introducing an 
off-diagonal correction to the Hamiltonian ^ 
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where are some complex constants of the order of unity. 
If \x\ < {hv/A) X {K/5K)^^-^/'^\ then each doublet in 
the spectrum of the Ising Hamiltonian ([g]) will split with 
the energy of the splitting 
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small compared to \5E\. 

Consider now the combined dynamics of electrons and 
spins at energies E <^ mui{hv/x, K). It is described by 
the Hamiltonian Hcg = Ho + Hj + A'H + U with 



U = hv 



Si'^\Xi){ii ■ T)'lp{xi) 
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where = [Puj{xi)K^]n. We retain only the x and y 
components of which cause electron backscattering, 



=Ci ■ {cos2kpXi,sm2kpXi,0) . 
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At temperatures T ^ e we may neglect the term AH in 
HcS- In that approximation, variables Si are constants 
of motion. For each configuration of Si we evaluate the 
backscattering current within the Born approximation in 
U. The total backscattering current is the Gibbs average 
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of such contributions. It leads to the correction to the 
balhstic conductance, 



5G 
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The term 2(^]^ • ^2)'? here comes from the interference 
between the electron waves reflected by the two local 
magnetic moments. The factor ■ ^2) experiences the 
conventional Fabry-Perot oscillations as a function of the 
Fermi momentum kp with the period 2'k/\x\. The factor 
77 = tanh[(5i<^cos(2fci;'a;)/T] also oscillates with the same 
period; at low temperatures, T <^ SE, it rapidly changes 
between —1 and 1 each time cos2fci?a; changes sign. 

At the lowest energy scale, T < e, one has to account 
for AH. It lifts the degeneracy of the ground state and 
therefore prevents an electron with energy less than e 
from elastic backscattering within the helical edge. In 
order to investigate the electron conduction in this low- 
energy regime, we evaluate the backscattered current in a 
steady non-equilibrium state induced by the source-drain 
voltage. The application of the Fermi Golden Rule yields 

/oo 
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+ [f++ip)f-+{p) - f+-{p)f-ip)] tanh — \ (12) 



where 



/a/j(p) ^ tanh + « ^^2r^^ ) ' o:,f3^±. (13) 



In the linear regime {V -> 0) Eq. (12 1 predicts a correc- 
tion to the ideal conductance 



SG = -(e? + ^2')^ (f , 0) - 2r;(|i . ^,)F (f , ^) (14) 
with 

2cosh^(A)cos(2zA) 
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[cosh (w) + cosh(2A)] 



(15) 



The second argument of F{w, z) is small in the entire 
region T <C SE/hP', so we may set z = 0. The u; 3> 1 
asymptote F(w, 0) ^ e^™ implies that the backscatter- 
ing correction is exponentially suppressed at T <C e. In 
the opposite limit w Q the function _F(w,0) — !■ 1 in 
agreement with Eq. (11) obtained in the AH = ap- 
proximation. Hereinafter we assume T ^ e.^^ 

Next, we generalize the above considerations to a sys- 
tem of A'^ > 2 magnetic impurities statistically uniformly 
distributed with average density n — N /L along the 
edge of length L. The corresponding effective Hamil- 
tonian which allows one considering scattering of elec- 
trons with energies E < hvn has the form T-LoS = 
Ho + 'Hisint' + AH + U, where AH and U are defined 



by Eqs. ([t]) and ^ (with extension of the summation to 
N), and 

hv X ^ ^i^j cos 2kp{xi — Xj) , 



Rising = > ■ 2^^lZZZ^Z±^Z^_ZAAs^s, , (16) 



with defined in Eq. ( 10 ) 



The conductance correction evaluated in the Born ap- 
proximation is 

SG = 2_/^*^J' cos2fci^(xi - Xj){s,Sj) . (17) 
Here the spin correlation function is 



Trexp j-Hising/T) s^sj 
Trexp(-'Hising/T) 



(18) 



In a given sample and at given temperature the con- 
ductance correction SG, Eq. (17), depends on the Fermi 



momentum kp through the oscillatory factors in Eq. (17) 
and in "Hismg, see Eq. (16 1. In experiment, this depen- 



dence take form of random fluctuations of the conduc- 
tance as the Fermi energy of electrons at the edge is 
changed. The statistical properties of such fluctuations 
are fully encoded in the cumulants 



Gm = lim — — lnexp(A • SG) 



(19) 



where the overline represents the statistical average. We 
now demonstrate that the statistical properties of con- 
ductance fluctuations at the magnetically contaminated 
helical edge are drastically different from those in the 
case of the conventional quenched disorder in a one- 
dimensional conductor. 

First, we recall the structure of conductance fluctua- 
tions caused by an ensemble of weak quenched scatter- 
ers in a usual one-dimensional conductor. At temper- 
atures such that the coherence length and the thermal 
length hv/T are both greater than the sample length 
the conductance is temperature- independent. The elec- 
tron reflection amplitude is, in the Born approximation, 
a linear superposition of A^ random complex numbers, 
^ikpxi ^ where Xi is the position of the ith impurity. Con- 
sequently, in the large A^ limit the conductance correc- 
tion obeys the Rayleigh distribution, which is essentially 
non-Gaussian. In particular, for the second and third 

3/2 

cumulants one has G3 = —^G^ ■ 

In contrast, the fluctuations of conductance of the he- 
lical edge, Eq. (17), exhibit strong temperature depen- 



dence in the whole range of validity of the model ( 16 1 



We note that the model possesses an intrinsic tempera- 
ture scale 



hv^n 
4tt 
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defining a crossover from the high-temperature regime, 
where spins thermally disordered to the "spin glass" 
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FIG. 1: MC data for the statistics of the mesoscopic con- 
ductance fluctuations at different temperatures and system 
sizes. In panel (a) shown is the average conductance gi as a 
function of temperature. The asymptotes given by Eq. (231 
are shown for comparison as solid black lines. Panels (b) 
shows the second cumulant 512 of the conductance distribu- 
tion. Panel (c) shows the dependence of the first three cu- 
mulants at T = 0.5 TsG on the system size. The log-log plot 
shows a good fit with the Qm ~ A*'™^^ scaling. Panel (d) 
shows the skewness of the distribution of conductance values 
as a function of temperature. 



regime where spins develop long-range correlations across 
the sample. For T ^ Tsg the model can be investigated 
by means of the virial expansion. 

We assume that the coupling constants have a Gaus- 
sian distribution with the average ^ and introduce the 
normalized cumulants gm such that 



Gr. 
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The normalized cumulants are functions of the dimen- 
sionless temperature r = T/Tgc and the number of im- 
purities N = nL. At T 3> Tsg the spin-spin correlation 
function in Eq. (17) is given by 



(siSj) = tanh 



hv^i^j cos 2kp{xi — Xj) 
Ati\x.i - Xj\T 



(22) 



Substituting Eq. (22) in Eq. (17) and averaging over the 



uniform distribution of magnetic impurities we find 
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The logarithm in Eq. (23) is due to the l/|xi — Xj| depen- 
dence in the large-distance expansion of Eq. (22). The 



upper cutoff for this dependence is provided by the sys- 
tem size, N/n, while the lower cutoff is defined by the 
distance x at which the exchange energy hv^'^/x equals 
the temperature. For the higher cumulants the virial ex- 
pansion yields 



gmiT,N) = Cr 



N 



m > 1. 



(24) 
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FIG. 2: Spin-spin correlation function at the opposite ends of 
the sample as a function of temperature for various systems 
sizes. Note that larger system sizes result in stronger cor- 
relations at given temperature. This effect is caused by the 
long-rangedness of the RKKY exchange resulting in the loga- 
rithmic renormalization of the spin-spin interaction constant 
with increasing system size. 



Here C„i are constants, in particular C2 = 0.20304 and 
C3 = 0.01242. Note, that unlike the quenched case the 
third cumulant is positive. At A^Tsq > T > Tsg the 
higher cumulants satisfy gf^-^ ^ 5™ therefore the distri- 
bution of conductance fluctuations is close to Gaussian. 

With decreasing temperature virial corrections to the 
spin-spin correlation function (18) become increasingly 
important. 



To explore this effect we calculate the first 
two terms in the virial expansion of the correlator 
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2T'^-n?\xi 



1 H — \a.Tn\xi — Xj \ + 



where . . . stand for the higher-order terms in 1/r. One 
can see that no matter how large the temperature t is, 
the virial expansion breaks down at sufficiently large dis- 
tances. For a given system size L one can define the 
crossover temperature such that 1 = r^^ln(nLT*). 



Below this temperature equations (23) and ([24| are not 
valid. The crossover temperature is a monotonically 
increasing function of the system size. In large systems 
there exists a parametric window 1 < r < where 
the long-range spin-glass correlations are absent, yet the 
virial expansion is invalid. 

In order to investigate g^ [t, N) at t < t» we per- 
formed Monte Carlo simulations for A^ — 25, 50, 100 and 
200. For each N we consider 10 random realizations of 
quenched impurity positions, assuming ^ to be the same 
for all impurities. We observe a considerable slowdown 
of the convergence of Metropolis algorithm for T < 5 Tsg 
caused by the onset of spin glass correlations. To over- 
come this difficulty we employ parallel tempering, which 
works efficiently down to T = 0.5 TgQ. Numerical results 
for 9i,2iT, N) are presented in Fig. [I] In the temperature 
window 5 < T < 10 the cumulants experience a sharp 
increase with decreasing temperature from g,„ ~ A^ (see 
Eq. ([24|) to 5,„ - TV^+i (see Fig.[l](c) ). For r < 5 one 



can observe the onset of correlations between moments 
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located near the opposite ends of the edge, Fig. [2j Note 
that at all temperatures the magnitude of the skewness 
of the conductance distribution is less than the universal 
value 1\pl predicted by the Rayleigh distribution. More- 
over, at both high and low temperatures away from the 
"spin glass" transition the skewness is suppressed indicat- 
ing a symmetric distribution of conductance, quite unlike 
the quenched case. 

To conclude, the purpose of this work is to reconcile the 
possibility of mesoscopic fluctuations in the conductance 
of a helical edge with the absence of coherent backscatter- 
ing in the presence of time-reversal symmetry (no exter- 
nal magnetic field applied). We find that scattering off 
an ensemble of large-spin (S > 1) magnetic impurities 
may open a temperature window in which the conduc- 



tance fluctuations are appreciable. The existence of such 
window is provided by a relatively strong effect of single- 
ion anisotropy which prevents easy flips of the impurity 
spins. It is further enhanced by the RKKY interaction 
between the spins. The latter interaction depends on the 
Fermi momentum of helical edge, bringing ergodicity in 
the conductance fluctuations as a function of the helical 
edge chemical potential. We elucidated the signatures of 
the described mechanism in the distribution function of 
conductance fluctuations. 
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Similar to 5G, one may consider the differential conduc- 
tance at r = 0. The differential conductance correction 
— 97bs / dV is exponentially small at e V <C e and oscillates 
with kp at a fixed eV > e. 



